Pinning effect and QPT-like behavior for two particles confined by a core-shell potential 
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We study the ground state entanglement, energy and fidelities of a two-electron system bounded by a core- 
shell potential, where the core width is varied continuously until it eventually vanishes. This simple system 
displays a rich and complex behavior: as the core width is varied, this system is characterized by two peculiar 
transitions where, for different reasons, it displays characteristics similar to a few-particle quantum phase tran- 
sition. The first occurrence corresponds to something akin to a second order quantum phase transition, while the 
second transition is marked by a discontinuity, with respect to the driving parameter, in the first derivatives of 
quantities like energy and entanglement. The study of this system allows to shed light on the sudden variation 
of entanglement and energy observed in Ref. [J. We also compare the core-shell system with a system where a 
core well is absent: this shows that, even when extremely narrow, the core well has a relevant 'pinning' effect. 
Interestingly, depending on the potential symmetry, the pinning of the wavefunction may either halve or double 
the system entanglement (with respect to the no-core-well system) when the ground state is already bounded 
to the outer (shell) well. In the process we discuss the system fidelity and show the usefulness of considering 
the particle density fidelity as opposed to the more commonly used - but much more difficult to access - wave- 
function fidelity. In particular we demonstrate that - for ground-states with nodeless spatial wavefunctions - the 
particle density fidelity is zero if and only if the wavefunction fidelity is zero. 



I. INTRODUCTION 

The realization of the importance of entanglement trig- 
gered a rethink in the way one can understand and quantify 
some quantum processes. Indeed, quantum information the- 
ory (QIT) has stemmed from the application of entanglement 
and the superposition principle to the processing and transmis- 
sion of data, 1 2] and it is now acknowledged that entanglement 
can play a central role in the description and understanding of 
quantum phase transitions (QPTs).[3-5] In QIT and QPTs it 
is important to determine how a quantum state changes under 
quantum operations or by varying external parameters. The fi- 
delity 12, 6] - extensively used in QIT to assess the 'closeness' 
of different quantum states - may naturally encompass the ef- 
fect of a driving parameter on a system, and, as such, it has 
been proposed as a key tool in understanding QPTs.|0-01 En- 
tanglement and fidelity can then provide a common language 
for QIT and QPTs.floWll The definition of a QPT has been 
widened by some authors to include changes in the quantum 
state of few-particle systems such as singlet-triplet transitions 
in a single quantum dot. lfl3ll Few -particle systems have also 
been used to characterize the predictive power of QPT indica- 
tors for a system undergoing a QPT in the thermodynamic al 
limit. (H 

In previous work, Ul it was shown that the transition from 
a core-shell to a double well potential induces a sudden varia- 
tion of both the entanglement and the energy of two electrons 
initially confined within the core well. This variation becomes 
sharper as the confining potential becomes harder, i.e. more 
similar to a rectangular-like potentials. This steep variation 
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was regarded as something potentially akin to a QPT but in 
the few-particle case. 

In order to understand this phenomenon, in this paper we 
will study systems related to Ref. [J and characterized by 
rectangular-like confining potential. We will focus on how 
ground-state entanglement, energy, and fidelities are affected 
by varying the potential core width and show that these simple 
systems encompass indeed a rich and complex behavior. The 
system we will mainly concentrate on is given by two elec- 
trons trapped within a core-shell potential, whose core reduces 
in width until it eventually disappears (see Fig.Q]). This may 
represent a (core-shell) quantum dot with an externally-driven 
confining potential: quantum dots are one of the most promis- 
ing hardware for the physical realization of QIT devices. Il5l- 
12311 hence, our findings may be of interest for QIT applica- 
tions. The system ground state is initially bound to the core 
well, but will become bound to the outer well (or shell) as 
the core width is reduced to zero and the outer well width 
increases. We will show that the corresponding sharp entan- 
glement variation is characterized by two very different transi- 
tions. The first presents elements akin to a second-order QPT 
and is associated with the transition of the ground state from 
the core to the outer well; the second is marked by a disconti- 
nuity in the energy and entanglement derivatives with respect 
to the driving parameter, and we demonstrate that it is due to 
the peculiarities of the confining potential. We will also ex- 
plicitly discuss the implications of these findings for the sys- 
tem described in Ref. Q]. 

Our analysis is important in the context of local sensitivity 
analysis. In particular, due to the pivotal role that the entangle- 
ment plays in severally quantum protocols (such as quantum 
algorithms !^ , quantum teleportation ll26ll and some quantum 
cryptography protocols 1I27I0 here we report on the sensitiv- 
ity of the entanglement with respect to small variations of the 
external parameter. ll28l l29ll characterize the region of the pa- 
rameter space over which the entanglement shows the steepest 



2 



variation and, consequently, ascertain the possibility of em- 
ploying the potential variations as entanglement 'switch'. In 
fact our calculations show that the presence of an inner core 
has a strong pinning effect on the entanglement even when the 
ground state is already bound to the outer well. Depending 
on the system geometry, it may in fact either halve the entan- 
glement (symmetric system) or double its value (asymmetric 
systems) when compared to the corresponding system without 
core well. This might potentially be exploited to induce sharp 
variations (switch) of the entanglement by modifying small 
regions of the confining potential. 

Finally, in the spirit of density-functional theory. ll30ll we 
will study whether the particle density can be used to track 
the system's ground state behavior via a particle-density fi- 
delity. We note that, from an experimental point of view, 
the density is a more accessible quantity than the full system 
wavefunction; our results show that, at least for the system at 
hand, the particle-density fidelity delivers information similar 
to the wavefunction fidelity. Importantly we will demonstrate 
that, for ground-states with nodeless spatial wavefunctions, 
the particle-density fidelity is zero if and only if the wavefunc- 
tion fidelity is zero. 



II. SYMMETRIC POTENTIAL, MODEL SYSTEMS 

We will first concentrate on systems with a symmetric con- 
fining potential (see Fig. [TJ. 

We consider three one-dimensional systems, each consist- 
ing of two interacting electrons bound by a confining potential 
and whose Hamiltonian in effective atomic units is 
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Here we set U (jgi,^) = 8{x\ — X2) to represent a con- 
tact Coulomb repulsion between the electrons. Vi(xj,R) 
are the confining potentials characterizing the three systems, 
i = DIW, OWO and DW, see below. 



A. System with a 'disappearing' inner well 
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Vdiw has a compact representation through the Heaviside 
step function, 



DIW 



(x;R) = V iw (x;R) + V ow (x;R), 



(6) 



where 



V iw (x; R) = V [9(x+(w-R)/2)6(-x + (w- R) /2)] 

(7) 

and V ow (x;R) = V %w (x;—R) describe the inner and the 
outer well, respectively. Eq. (|6]l is equivalent to Eqs. (|4]i and 
if we assign 9(0) = 0. This is consistent with considering 
the Heaviside step function 6{x) as, for example, the limit (in 
a distribution sense rf3lll ) forp — ► 00 of 
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where p and m are positive integers. With p ~ 10 and 
m ~ 20, we get a smooth, 'softer' version of Vdiw- As 
p — > 00 arguments similar to the ones developed in Ref. |T| 
seem to suggest a discontinuity in the entanglement entropy 
and energy derivatives (and hence something reminiscent of a 
QPT in the few-particle regime). The chosen parametrization 
for the potential will help us to better understand this limit. 



B. Benchmark system 

The confining potential of the 'outer well only' (OWO) sys- 
tem is given by Vqwo = V ow , see inset of Fig.Q] We use this 
system as a benchmark. 



The potential of the 'disappearing' inner well (DIW) sys- 
tem, Vdiw(x', R), is characterized by an inner (core) and an 
outer shell well, see Fig.Q] As the parameter R increases, the 
inner well width, W lw , becomes narrower and the outer well 
width, W ow , increases as 



w 
0. 



R for R < w 
for R > w 



W™{R) = 
W ow {R)=w + R. 
Taking Vo as the depth of the outer well, we can write 



(2) 
(3) 



C. Core-shell to double well system 

This is the rectangular-like limit of the system considered 
in Ref. [lj. As the driving parameter changes, this potential is 
modified from a core-shell to a double-well potential. The ex- 
plicit expression of this potential in the rectangular-like limit 
can be written as 

Vdw(x;R) = V o [0(-x + (w - R)/2)6{x + (3w - R)/2) 

+6{-x + {3w - R)/2)6(x +{w- R)/2)] .(9) 
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Here the transformations Row = 2iu — R and d = w/2 give 
the control parameter and the inter-well distance as used in 
Ref.Ql respectively. 

For the subsequent calculations, unless otherwise stated, we 
use w = 5 ao, where ao is the Bohr radius, and Vq = —10 
Hartree. 
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FIG. 1: Potential Vdiw versus x for R = 4ao. Inset: same as main 
panel but for R — 5 ao, for which DIW and OWO systems coincide. 



HI. RESULTS FOR ENTANGLEMENT AND ENERGY 
(DIW AND OWO SYSTEMS) 

To calculate the ground-state properties, we directly diago- 
nalize the Hamiltonian Eq. (Q~|i, by writing its eigenfunctions 
$t as a linear combination of single-particle basis functions 
and truncating the corresponding expansion as 



M M 

jl = l j2 = l 



2 - k r] jl (xi;uj)T] j2 (x2;uj), (10) 



where r\. (x; u>) are the eigenfunctions of the one-dimensional 
harmonic oscillator with angular frequency u. A single- 
particle basis size of M = 50 with cj = 2 ensures good con- 
vergence of the results at any R. 

We calculate the particle-particle spatial entanglement|32] 
and the ground-state energy of the system for 4 ao < R < 
8 ao. For the DIW system R = 4 ao corresponds to a core- 
shell structure with the two electrons confined in the inner 
well, while for R > 5 ao we have Vdiw — Vqwo- 



A. Energy 

First we consider the ground state energy Eq of the DIW 
system (solid line in Fig. |2}\) against the benchmark (OWO, 
dashed line). 

As R becomes larger, the inner well narrows and the en- 
ergy of the two-electron state increases, until the electrons are 
eventually 'forced' into the outer well. The ground state en- 
ergy leaves the inner well at R = R c = 4.96 ao. This cor- 
responds to an inner to outer well ratio of 0.0039. Hereafter, 



we will refer to the parameter region around R c as the 'mi- 
gration region' : for these values of the driving parameter the 
system wavefunction is the most sensitive to driving param- 
eter changes. Here the electron wavefunction 'expands' into 
the outer well and, as a consequence of this, the system shows 
the most interesting behavior. This region of high sensitivity 
is relatively narrow and in fact for R ^ 5 ao the ground state 
energy becomes a very slowly-varying, decreasing function of 
R. 

The first derivative of the ground state energy with respect 
to the driving parameter, dEo/dR, displays a discontinuity 
at R = 5 ao, but it is smooth elsewhere (see Fig. |2j}). This 
discontinuity is found in the first derivatives with respect to 
R of all the quantities we consider. dEo/dR has a maxi- 
mum at R = 4.7 ao. From Fig. [2)3 (inset and main panel) 
we see that at first the shrinking of the inner well increases 
the ground state energy with an increasing "speed". However, 
in the migration region the change in the ground-state energy 
rapidly slows down: in this region the wavefunction is starting 
to spread into the larger outer well, hence moving towards a 
regime where Eq is almost constant with R. 

The second derivative of Eq with respect to R displays a 
marked minimum at R = 4.90 ao and an infinite discontinuity 
at R = 5 ao, see Fig.[2f . 

The behaviors of the Coulomb energy (U), and of the ki- 
netic energy (T) are plotted in the upper panel of Fig. [3] 
where (...) indicates the ground-state expectation value. For 
the DIW potential, both display a maximum located at R = 
4.47 ao (corresponding to an inner to outer well ratio of 
0.058). The ratio between the Coulomb and the kinetic in- 
teractions, Fig. [3)3, provides an unambiguous signature of the 
migration point R c , whereas no particular structure emerges 
from the visual inspection of both Coulomb and kinetic en- 
ergy separately, Fig.[3|\. 



B. Entanglement 

We calculate the spatial entanglement j32| using the von 
Neumann entropy S and the linear entropy L, 

S = -Trp red log 2 p red , (11) 
L = Tr(p red - p 2 red ) = 1 - Trp 2 red , (12) 

where p red = Tta \^f) (^ \ is the reduced density matrix found 
by tracing out the spatial degrees of freedom of one of the 
two particles (subsystem 'A' ) and ^ is the ground-state. We 
consider also the position space-information entropy S n , 



n(x) hin(x)dx, 



(13) 



where n(x) is the system particle density. 

For a pure bipartite state the von Neumann entropy S is the 
unique function that satisfies all the entanglement measure- 
ment conditions. ll33l [34tl while the linear entropy L is com- 
putationally convenient and quantifies the entanglement in the 
sense that it gives an indication of the number and spread of 




FIG. 2: Ground state energy Eo (panel A), first and second deriva- 
tive of -Bo with respect to R (panel B and C, respectively) for the 
DIW (solid line) and OWO (dashed line) potentials as functions of 
the driving parameter R. In all the three panels the inset zooms on the 
'migration region' with R c — 4.96 <xo indicated by a vertical dotted 
line. 



terms in the Schmidt decomposition of the state. The position- 
space information entropy S n can be considered as an approx- 
imation to S when off diagonal terms are neglected 0211 and 
is written in terms of the particle density, so it could be more 
easily and directly accessed by experiments. 

In Fig. 03, the linear, von Neumann, and position-space 
information entropy are plotted as a function of R for the 
DIW system. L and S n have been rescaled so that they have 
the same value of S at R = 8 do- All quantities show the 
same qualitative behavior, L and S rescaling almost perfectly 
onto each other. In particular, all quantities show a non- 
differentiable point at R = 5 ao and present a minimum lo- 
cated in the same R region. However, the minimum of S n 
(R = 4.45 ao) is nearer to the maximum of (U) than the 
minima of the other two entropies (R = 4.51 ao for L and 
R = 4.52 ao for S), and is more pronounced. 

By the Hohenberg-Kohn theorem, Ollll the ground state par- 




FIG. 3: Panel A: Coulomb energy, {[/), and kinetic energy, (T), for 
the DIW and OWO potentials as a function of the driving parameter 
R. Inset: as main panel, but in the neighborhood of the migration 
point R c , marked by a vertical dotted line. Panel B: Ratio between 
the Coulomb interaction energy and the kinetic energy, (U)/{T), 
versus R for both the DIW and OWO potential. Inset: details of the 
'migration region' with the vertical dotted line indicating the point 
R — R c - Panel C: The von Neumann (S, dotted line), and rescaled 
linear (L r , solid line) and space-information (S n>r , dashed line) en- 
tropies as functions of R for the DIW system. The rescaling was 
chosen in such a way that L r and S„ tT are equal to S at R — 8 ao. 
This results in L r — 3.041/ and S„ r = 0.15S n . 



ticle density uniquely determines all the ground-state proper- 
ties of the system, so in principle the ground-state entangle- 
ment for this system could be written as a functional of the 
density; the overall similarity between S n - explicitly written 
as a functional of the density - and the two entanglement mea- 
sures S and L reinforces the idea that pertinent information 
can be extracted from the electron density. As for the DIW 
system L can be rescaled very well onto S, we will continue 
using the computationally convenient linear entropy L. 

The linear entropy of the DIW, Ldiw> an d of the OWO sys- 
tem are compared in Fig. |4j\. Ldiw displays three regions. 
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FIG. 4: Panel A: Linear entropy (L) as a function of R for the DIW 
and OWO potentials. The dots indicates the value of L diw at R = 
R c . Panels B and C: first and second derivatives with respect to R of 
the the linear entropy as a function of R. In all three panels, the inset 
represents the respective function around the migration point R c , the 
latter being highlighted by the vertical dotted line. 



of (U) corresponds here to a minimum of the entanglement as 
these extrema occur when the electrons are most confined and 
hence in an almost factorized state. [[ll] The decrease of (U) 
is a signature of the wavefunction spilling into the outer well 
and, consequently, of an increasing influence of Coulomb cor- 
relations in shaping the wavefunction with a corresponding 
increase of the entanglement. In the migration region (with 
(U) (R c ) being approximately 7% of its maximum value), the 
wavefunction density is substantially spread within the outer 
well, and small variations of R produce large changes in the 
entanglement. 

We note that in Ref.[35t a system similar to DIW, with an in- 
ner well shrinking in width but never disappearing, was stud- 
ied. In this a case no discontinuity in any derivative of the 
relevant quantities were found. 



C. Comparison between the DIW and the OWO potentials 

In Figs. EH andgK E and L are plotted for both the OWO 
and DIW potential. At R = R c the many-body ground-state 
is bounded to the outer well, and in particular E® IW (R c ) ss 
0.99Eg wo (R c ). In contrast, L DIW (R C ), which is marked 
by a dot in Fig.Hl is approximately half of the corresponding 
entanglement value in the OWO case. We underline that here 
the inner well has a finite depth but the inner to outer well ratio 
is only 0.0039, so, from a geometric point of view, the inner 
well should be negligible. However in the region R c < R < 5 
the behavior of (U), with (U) DIW > {U) owo , suggests that 
the electrons remain strongly pinned to the inner well region 
even though the width of the latter is basically negligible. As 
a consequence a very narrow inner well is able to modify the 
distribution of the electrons in such a way that their entangle- 
ment is highly and non-linearly reduced. This 'pinning prop- 
erty' of the entanglement might open possibilities of rapidly 
and efficiently modifying the entanglement in a nanostructure 
system. 



IV. GROUND-STATE WAVEFUNCTION AND 
PARTICLE-DENSITY BEHAVIOR 



The first is characterized by a slow variation in entropy with 
a shallow minimum at R ', = 4.51 an. In the second region 
(4.9 ao < R < 5 a) the entropy increases very rapidly, and fi- 
nally for R > 5o,q the entropy increases linearly with R. The 
first derivative of Lbiw (Fig. |4jB) presents a shoulder-like 
structure connecting the first and second regions; then, after 
the boost in the rate of change of the entropy, the derivative has 
a finite discontinuity at R = 5ao. The second derivative of 
Ldiw presents two maxima at R = 4.89 ao and R = 4.99 ao, 
and a minimum at R = 4.92 ao. It has an infinite discontinuity 
at R = 5 ao 

The rate of change of the entropy shown in Fig. [4] is the 
result of the competing effects of the confinement strength 
and of the Coulomb repulsion. However, since it is the ra- 
tio between these two factors which governs the response of 
the system to a variation of the driving parameter, a maximum 



In Fig.[5]the high sensitivity of the system wavefunction to 
small changes of the driving parameter in the migration re- 
gion is explicitly demonstrated. The figure in fact shows the 
wavefunction contour plots for R = 4.5 ao (~ maximum of 
(U) DIW , panel A), R = R c ('migration' point, panel B) and 
R = 5 ao (value at which the inner well disappears, panel 
C). The wavefunction becomes more and more confined until 
R » 4.5 ao, for which it displays a single maximum (panel 
A). A further reduction of the inner well width induces the 
wavefunction to leak into the outer well (compare scales on 
axis of panels A and B). Around R w R c the wavefunction 
starts to separate into two lobes, but remains largest close to 
the inner well (pinning effect). As R increases beyond R c , the 
shape of the wavefunction displays two well-defined lobes, 
reflecting the effect of the electron-electron repulsion com- 
bined with the diminished confinement strength (panel C). 
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The wavefunction width and height though remain roughly 
constant, compare panels B and C. We note that the wave- 
function shape appears to change "smoothly" as R increases 
and, in particular, no detectable change in the geometry of 
the wavefunction seems to take place at R = 5 <xrj, where the 
non-differentiable points of the entropy and energies are both 
located. 
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FIG. 5: (color online) Contour plot of the wavefunction against the 
particles' positions x\ and X2 for the DIW potential at R = 4.5 ao 
(~ maximum of (U) DIW and minimum of L, panel A), R — R c — 
4.96 ao (panel B), and R = 5 ao (point at which Wo" IW = 0, panel 
C). 



As already mentioned, the particle density n(x;R) should 
uniquely capture the system ground-state behavior, so we will 
now check if the density shape is more susceptible than the 
wavefunction to the shrinking and disappearance of the inner 
well. In Fig. [6] the density is plotted for various values of 
R. We see that, as R increases, the height of the central (and 
only) peak diminishes. For R s=s R c the density develops two 
shoulders and at R = 5 ao the central peak disappears and is 
replaced by two peaks which are symmetric around the ori- 
gin. This can more clearly be seen in the inset. At least for 
the system at hand, the pinning from the inner well has a more 
clear-cut effect on the shape of the density than on the shape 
of the wavefunction, as in particular it determines the pres- 
ence or absence of a central peak for the particle density. At 



difference with the wavefunction, the change in the number of 
peaks of the particle density associated to the disappearance of 
the central maximum can then be associated with the discon- 
tinuity in the derivatives of energy and entanglement caused 
by the disappearance of the inner well. 




FIG. 6: Density n(x; R) for the DIW potential plotted against the 
position x for four different values of R (as labeled). Inset: zoom of 
main panel for R « 5. 



V. FIDELITY OF THE GROUND-STATE WAVEFUNCTION 

The fidelity between two states quantifies their similarity 
and as such has been extensively used in quantum information 
theory. [|2j] More recently the fidelity has been introduced as a 
method for the characterization of OPTs:lr7[l8[ l36l] a signature 
of QPT is an abrupt change in the wavefunction. 0711 and this 
suggests the evaluation of the fidelity between states across 
the critical point as a good choice for the identification of a 
QPT. Here we will use this method to better understand the 
system behavior in the migration region. 

For our system the ground-state fidelity is given by 

F(R 1 ,R 2 ) = | {if>(xi,X2-,Ri) \ip{x u x 2 -R 2 ))\- (14) 

following Eq. (TTOb we then calculate it as F(Ri,R 2 ) = 

^2ji,j2=i a h,h,n(Ri)aji,j 2 ,n(R2) ■ F(R 1 ,R 2 ) can be in- 
terpreted in two complementary wavs. l38ll and according to 
this we will consider two different sets of values for R\ and 
R 2 . 

In quantum information theory, the fidelity can be seen as a 
generalization of a measure of similarity between two classi- 
cal probability distributions. [2] Let us take Ri = Rq, where 
ip(Ro) is the reference state, then the fidelity is the overlap 
between this initial state and the wavefunction ip(R) calcu- 
lated as R 2 = R varies in the parameter space. The fidelity 
F(Rq, R) clearly depends on the choice of the reference state. 
The fact that the minimum of the entanglement corresponds to 
a quasi-product state (see Fig.|H panel A), which evolves to- 
wards a highly entangled state as R increases, suggests as a 
natural choice Rq = 4.52 ao, corresponding to the minimum 
of the linear entropy. 

Alternatively, the fidelity can be seen as a geometrical ob- 
ject connected to the Fubini-Study distance between quan- 
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turn states, H] where the square distance between infinitesi- 
mally close states can be approximated as ds^, s ss 2(1 — F). 
In this case the fidelity is calculated between two wavefunc- 
tions depending on infinitesimally different parameters, i/j(R) 
and i/j(R + SR). At the critical point, where there is an 
abrupt change in ip, this function has a minimum and pos- 
sibly a discontinuity. In Fig. |7J the fidelities F(Rq,R) and 



functions and N particles can be approximated as 

F(R, R+SR) « 1-^/ ( d ^ dR X ^ R) y dxi . . . dx N , 

(15) 

shows a discontinuity at R = 5 ao, in accordance with the 
discontinuity found in the derivatives of all the quantities dis- 
cussed so far. 




4.5 4.7 4.9 5.1 5.3 5.5 5.7 5.9 6 

R (a.) 




4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8 6 

R C„) 



FIG. 7: Panel A: fidelity F(R , R) vs R for i? = 4.52 a . Inset: 
first derivative of F(R , R) vs i?. Panel B: fidelity F (R, R + 5R) 
vs R with SR = 0.002. Inset: as for the main panel, but zooming on 
the region of the fidelity minimum. The plots refer to the DIW and 
OWO potential (as labeled); the vertical dotted line corresponds to 
Re 



F(R, R + 5R) are plotted as a function of R (panel A and 
B, respectively). F(Rq, R) displays three distinct regimes, in 
accordance with the behavior of all the other quantities stud- 
ied so far. In particular, we see that for 4.8 clq < R < 5 clq 
we have a dramatic decrease of the fidelity: the wavefunc- 
tion is rapidly changing from a quasi-product state towards 
a triplet-like entangled state yj] (compare Fig. |3J panels A 
and C). The derivative dF(Ro,R)/dR presents a minimum 
at R « 4.92 ao, near the migration point R c . For R > 5 ao 
the fidelity is almost constant and drastically reduced, with 
F(Ro,R = 5) ~ 0.19: in this region the wavefunction 
is nearly orthogonal to the reference state. We note that 
F(Ro, R) is not differentiable at R = 5 ao. 

The behavior of F(R, R + 5R) (Fig. |7fc) shows that the 
most significant changes in the wavefunction are confined to 
the migration region, around a marked minimum at R = 4.93, 
again very close to R c . F(R, R + SR), which for real wave- 



VI. FIDELITY OF THE PARTICLE DENSITY 

For a 1-particle system with control parameter R we have 
that the particle density is n(x; R) = \ip(x\R)\ 2 so, in this 
case, the fidelity may be written in terms of the density as 
F(R 1 ,R 2 ) = J y /n(x;R 1 )n(x;R 2 )dx. 

We may generalize this to a 'density fidelity' by using the 
density arising from A^-particle systems, 

n(x,R) = N J \ip(x, x 2 ■ ■ ■ , x n ; R)\ 2 dx2 ■ ■ ■ dx n . (16) 
and defining the 'density fidelity' as 

F n (R 1 ,R 2 ) = jj J y/n(x;R 1 )n(x;R 2 )dx. (17) 

F n (Ri, R 2 ) has the properties expected from a fidelity, that 
is < F n (Ri, R 2 ) < 1 and it measures the overlap between 
particle densities as the driving parameter R is varied. We will 
also demonstrate that F n (R±, R 2 ) vanishes if and only if the 
corresponding wavefunction fidelity F(R±, R 2 ) vanishes. 

We note that a density fidelity has been proposed for lat- 
tice systems and linked with QPTs in Ref. [33- We initially 
calculate the density fidelity with respect to i?o- F n (Ro,R) 
shows a non-differentiable point at R = 5 ao corresponding 
to the disappearance of the inner well (Fig. |8}\); its deriva- 
tive in respect to R is plotted in the inset. We note the simi- 
larity between the behavior of F n (Ro, R) and F(Ro, R) and 
between their derivatives, the main difference being that the 
residual fidelity for R > 5ao is larger for the density than for 
the wavefunction. 

We show F n (R, R + SR), where SR = 0.002, in Fig. 1 
(panel B). As for the corresponding wavefunction fidelity, 
here the discontinuity at R = 5 ao appears directly in the 
'density fidelity'. Again the behavior of F n (R, R + SR) and 
F(R, R + SR) are very similar (compare Fig. [8^ with the in- 
set of Fig. |7jJ), with the density preserving a slightly higher 
fidelity at its minimum, which occurs at R = 4.92 ao 

In the DIW system, viewing the density seems to more 
clearly and readily display the fast changes in the ground state 
properties corresponding to the discontinuity in the derivatives 
of Eq and L than viewing the wavefunction (see comments to 
Figs. [5] and [6]). This may be due to the lesser formal com- 
plexity of the density, which is always a function of a sin- 
gle position vector - x in the present case - as opposed to 
the complex many-body wavefunction, a function of N posi- 
tion vectors whose parameter space is clearly more difficult 
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FIG. 8: Panel A: Density fidelity F n (Ro, R), Ro = 4.52 a , for 
the DIW system, plotted against the parameter R. Inset: Derivative 
of F n (Ro, R) with respect to R plotted against R (Ro — 4.52 ao). 
Panel B: F n (R, R + SR) with SR = 0.002, plotted against the pa- 
rameter R for the DIW system. 



to analyse and visualize. In addition the particle density fi- 
delity is able to predict all the other notable features of the 
wavefunction fidelity, such as the minimum occurring around 
R « R c . As noted, minima in the fidelity F(R, R + 5R) are 
associated to abrupt changes in the wavefunction and may sig- 
nal the occurrence of a QPT, so, in accordance with Ref. [39l 
our results suggests that the density fidelity may be used as 
an alternative to the wavefunction fidelity to understand brisk 
changes in the ground state and hence to study QPTs. This 
is in line with the Hohenberg-Kohn theorem which in its sim- 
plest form shows that for non-degenerate ground-states, the 
density uniquely determines the many-body wavefunction and 
so all the properties of the system. [13011 We point out that the 
particle density is a much easier quantity to calculate (and to 
experimentally access) than the full many -body wavefunction. 
As such the use of the fidelity density might become of great 
help in understanding phenomena such as QPTs. Similarly its 
characteristics as highlighted above suggest that it could be a 
useful tool for local sensitivity analysis. 



A. One-to-one correspondence between vanishing of ground 
state particle density fidelity and wavefunction fidelity 

We will now demonstrate the important property that, for 
systems with finite external potentials and ground-states with 
nodeless spatial wavefunctions, the density fidelity is zero if 
and only if the ground-state spatial wavefunction fidelity is 
zero. 

The nodeless spatial ground-state wavefunction of a time- 
independent Hamiltonian may always be taken to be real and 



positive, so for any two such A^-particle ground-state wave- 
functions -01 and ip2 we may define the real positive function 



. dx 



N- 



(18) 



For any fixed x this defines an inner product, as positive def- 
initeness is satisfied by f\ : i(x) = n(x)/N > for finite ex- 
ternal potentials. The Cauchy-Schwarz inequality can then be 
written as 



/ ipiip2dx 2 ■ ■ - dXN 
< (j \ipi\ 2 dx2 ■ ■ ■ dxN J \?p2\ 2 dx 2 ■ ■ ■ dxN 
integrating both sides with respect to x leads to 

F < F„ 



(19) 



(20) 



and so if F n tends to zero, so must F. 

For a general wavefunction, a fidelity of zero does not 
imply a density fidelity of zero as for example two excited 
state wavefunctions may both be non-zero in some finite re- 
gion of space but still be orthogonal. In addition, if we 
compare wavefunctions arising from different forms of inter- 
particle interactions, say attractive and repulsive, then, again, 
the density cannot always discriminate between orthogonal 
wavefunctions. This can be explicitly seen by considering 
the limiting case of infinite inter-particle attraction or repul- 
sion. Let us consider two particles in one dimension: in 
the case of infinite attraction their wavefunction will sat- 
isfy ipAixi, x 2 ) = if X\ 7^ x 2 , while for infinite repul- 
sion we have ij}R{ x i> x 2) = if X\ = x 2 , otherwise both 
%p > 0. Clearly we obtain J ipA'4'Rdxidx 2 = 0. Let us 
now consider the single particle densities. In general it will 
be ur(x) = f \ipR(x, xi)\ 2 dxi > 0. To ensure normal- 
ization, iPa{xi,x 2 ) 2 = <f)i(xi) 2 5(xi — x 2 ), with (f)\(x\) it- 
self normalized, giving a corresponding single particle density 
71,4(2;) = |^i (x)| 2 > 0. It follows that the related density fi- 
delity J yj njx{x)riA(x)dx is different from zero. 

However, if we assume the requirements needed for stan- 
dard DFT, i.e. ground state and same inter particle interaction, 
then we can argue that the density fidelity can detect orthogo- 
nal, nodeless ground-state wavefunctions. The lack of nodes 
in the ground states means that we can choose a phase so that 
both our wavefunctions are never negative. Here a fidelity of 
zero corresponds to the hypothetical situation when the wave- 
functions do not overlap at all. When the inter particle inter- 
action is fixed this lack of overlap arises because the wave- 
functions are spatially distinct, and so the densities will not 
overlap. Hence for ground-states with nodeless spatial wave- 
functions, the density fidelity is zero if and only if the spatial 
wavefunction fidelity is zero. 



Vfl. QPT-LfKE TRANSITION (SYMMETRIC SYSTEMS) 

As pointed out previously, a minimum in F(R, R + 8R) 
may highlight a QPT and certainly witnesses a rapid change in 
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the wavefunction. In our case the, minimum in F(R, R + 6R) 
observed in Fig. [7] corresponds to the transition between two 
separate sets of ground states; the first set bounded by the in- 
ner and the second set bounded by the outer well. Fig. |7j\ 
shows that this transition is between states that are almost or- 
thogonal. As the width of the inner well is reduced, the en- 
ergy gap between these set of states reduces: this transition 
has some of the characteristics of a second-order QPT. 

This is apparent when looking at the ground state energy 
derivatives: d 2 E n /dR 2 presents in this region a marked min- 
imum, which in turn corresponds to an inflection point in the 
energy first derivative. If this were a full-fledged QPT, this 
inflection point would have a vertical tangent, and hence the 
minimum in d 2 Eq / dR 2 would become a divergency. 

As discussed in Refsl40landl4ll a second-order QPT should 
be signaled by a corresponding structure in the first derivative 
of the entanglement. The first derivative of the entanglement 
entropy presents indeed a structure (a shoulder) whose width 
can be defined by the first maximum-minimum structure in 
d 2 L/dR 2 , i.e4.89a < R < 4.92 a (see Figg] panels B and 
C): this shoulder indeed frames the region of the minimum 
of d 2 Eo/dR 2 . The bulk of the wavefunction change should 
occur in the region of the minimum of F(R, R + SR): in 
Fig.|9]we then present the wavefunction at R = 4.91 cio (panel 
A) and R = R c = 4.96 ao (panel B). The plots confirm a 
quite substantial change in the wavefunction, which smears 
over the upper well as R increases, changing from a single, 
pointed peak towards a two-lobe geometry. 

As for the case of a finite-size system which would undergo 
a QPT in the thermodynamic limit, [42] the transition we ob- 
serve in the wavefunction occurs over a (small) parameter re- 
gion and slightly away from the expected 'critical' value of 
the driving parameter, i.e., for R < R c . 



DIW: R=4.91 



ORIGIN OF THE DISCONTINUITIES OBSERVED 

AT R — w 




X 2 ( a o) 4 



B)DIW: R = 4.96 «, 




FIG. 9: Ground-state wavefunction plotted against the particles' 
positions x\ and xi for the DIW potential at R — 4.91 ao , panel A, 
and R = R c = 4.96 ao, panel B. 



VIII. 



It was demonstrated in Ref. |4Cj |41| that a discontinuity in 
the first (second) derivative of the ground-state energy with 
respect to the driving parameter - a signal of a QPT - may 
correspond to a discontinuity in the (derivative of the) ground- 
state entanglement. 

What we observe in the present case at R = w is instead a 
discontinuity in the same order derivatives of the ground-state 
energy and entanglement. Moreover, all the other quantities 
under study, such as (U) or S n , are non-differentiable at the 
same point. A similar situation was speculated for the limit 
p — >• oo in Ref. [J Here we would like to clarify the origin of 
the discontinuities we observe. 

First of all we can extract from the fidelities important infor- 
mation on the ground state wavefunction behavior at R ', = w: 
the continuity of F(Rq, R) shows that the ground-state wave- 
function is continuous, on average, at R = w (see Fig. |7J\). 
However the discontinuity of 8F(Rq, R)/dR\n =w indicates 
a discontinuous derivative for the wave function at the same 
point (inset of Fig.[7J\). The discontinuity of dip(x; R)/dR at 
R = w is confirmed by the discontinuity of F(R, R + 5R) at 
the same point, see Eq. (TT3T >. 

To understand the above picture we consider the Hamilto- 
nian for a general potential 



H{R)=H +J2v(x j ,R), 



(21) 



where Hq = T + U is independent from the driving param- 
eter R, T is the kinetic energy, and U is the electron-electron 
interaction. From the Hellmann-Feynman theorem we have 



dE{R) 
dR 



dV{xj]R) 



dR 



i/>(x;R) 



(22) 



where x = (xi, . . .xjv) represents the coordinates of the N 
particles. Consequently, if dV(xj\R)/dR is discontinuous 
with respect to the parameter R, then this discontinuity could 
propagate to dE/dR. 

We note that in the case of a first order QPT, the disconti- 
nuity in dE/dR should arise from the wavefunction, and not 
from the potential. In the present case, while the fidelity indi- 
cates a continuous wavefunctions, not only dE/dR, but also 
all the other quantities used as indicators for a QPT, present a 
point of non-analyticity at R = w = 5 ao. It is hence neces- 
sary to understand how a discontinuity in the potential would 
affect the other quantities of interest, and in particular if, in 
contrast with the situation in Refs. |4l], E(| it would produce 
discontinuities in the first derivative of both energy and entan- 
glement entropies. 

By considering the time-independent Schrodinger equation 
associated to Eq. d2Tb we can write 



10 



TiP{x;R) UiJj(x;R) 
ip(x;R) ip(x; R) 



E(R) = J2v(x f ,R). (23) 

3=1 



Eq. (1231 is well-defined for a nodeless ground state wave- 
function and can be seen as a family of equations labeled by 
the continuous parameter R. Eq. d23l shows that if the poten- 
tial is discontinuous only at a set of points of measure zero 
then, at most, it may only directly cause ip and/or Tip to be 
discontinuous on that same set of points. Hence, for finite 
discontinuities, these discontinuities will not propagate to any 
integrated quantities such as expectation values. We then con- 
tinue by assuming that ip and Tip are, at worst, discontinuous 
over a set of points of measure zero. We differentiate Eq. 
with respect to R, and use Eq. ( 1221 to obtain 



-T 



dip{x;R) [Tip(x; R)} dip(x; R) 



ip(x;R) dR 
dV(xf,R) 



E 



dR 



ip 2 (x;R) dR 
dV{x j ;R) 



iP(x;R) 



dR 



1>(x;R) 



(24) 



Finite discontinuities in the potential may mean that its deriva- 
tive with respect to R will comprise delta functions, and 
hence that, unless accidental cancellations occur, these dis- 
continuities will propagate to (dV{xj]R)/dR). Let us as- 
sume that they are such that (dV(xj; R)/dR) is discontinu- 
ous at R = R. Then, on the right-hand side of Eq. d24l we 
have two discontinuous functions with respect to R, but as the 
second term does not depend on x, the right-hand side is ac- 
tually discontinuous at (x; R) for all or almost all values of 
x since no accidental cancellations can hold for all x. This 
means that the left hand side of Eq. ( 1241 will present the same 
discontinuities. As ip(x;R) and Tip are at least continuous 
almost everywhere with respect to x at (x; R), this implies 
that dip(x; R)/dR has indeed to be discontinuous at (x; R) 
for all or almost all x and hence the first derivative in respect 
to R of any functional of ip will be in general discontinuous at 
R = R. This is exactly what we observe. 

In Appendix Al we illustrate these points by explicitly 
analyzing the effect of the finite discontinuity at the point 
(x = 0; R = w) in Vdiw- 

In the next section we will instead consider a counter- 
example for which, due to an accidental cancellation, 
(dV(xj; R)/dR) - and hence all first derivatives in respect 
to R - remains continuous even in the presence of discontinu- 
ities in V(xj ; R) similar to the ones of the DIW potential. 



IX. CORE-SHELL TO DOUBLE WELL POTENTIAL 

In Ref. Q] it was speculated that, in the rectangular-like po- 
tential limit, the observed sharp transitions in energies and 
entanglement would display non-analyticities as the potential 
changes from a core-shell structure to a double well potential. 



The behavior of the entanglement and its derivatives in this 
limit is shown in Fig.fTOl 
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FIG. 10: Upper panel: Entanglement entropy L plotted against R for 
the potential of Eq. ((9} and w = 5ao. Inset: same as main panel 
but plotted in respect to Rdw for an easier comparison with Ref.Q]. 
Middle panel and lower panel: dL/dR and d 2 L/dR 2 , respectively, 
vs R for the potential of Eq. ^ and for w = 5 do. The insets present 
the transition region and the vertical dotted line the migration point 
R c = 4.98 a for the potential Eq. 10. 

For the DW potential Eq. ©, the plot of F(R 0l R) 
(Fig. [TTB. with Rq = 4.76ao the minimum of L for this 
system) confirms that the reference state, almost factorized 
and bounded to the inner well, is practically orthogonal to 
the triplet-like ground state reached after the transition to the 
double well potential.^ However a QPT-like transition oc- 
curs only when the ground state bounded to the inner well 
migrates to the outer well, see the shoulder in dL / dR and the 
minimum of d 2 Eo/dR 2 in Figs. [TOB andfTTK. respectively. 
The subsequent transition to a double-well potential merely 
further isolates the two lobes of the wavefunction from each 
other. The sharp increase in entanglement which corresponds 
to this further change, does not then signal any further QPT- 
like point, as confirmed by the absence of additional structures 
in d 2 E /dR 2 and F(R,R + 5R) (see Fig.QjK and Fig.Qjt)- 
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FIG. 11: Panel A, B and C: d 2 E /dR 2 , F(R , R) and F(R, R + 
SR) for the DW potential with respect to R. The insets in panel 
A and C present the details of the transition region; the inset in B 
depicts the first derivative of F(Ro, R) with respect to R and with 
Ro = 4.76a . 

Interestingly, in contrast with the DIW case, the disconti- 
nuity at R = w of (dV(xj] R)/dR) is accidentally removed 
in this system. As a result the derivatives with respect to R of 
the various quantities of interest, as well as F(R, R + SR), do 
not present discontinuities (see Figs. PTOl and [TTli. Details are 
given in Appendix A2. 



X. ASYMMETRIC POTENTIAL 

We will now consider the possibility for the position of the 
inner well, with respect to the symmetry axis of the outer 
well, to be shifted toward left (see Fig. [12] upper panels) and 
we will refer to these asymmetric potentials as Asymmetric 
Disappearing inner Well (ADW) potentials. In Fig [T2b, the 
inner well position is left-shifted by a quantity w a = 2 ao, 
where w a is the distance between the symmetry axes of the 
inner and outer well. In panels PT2C and [12)3, w a = 1.5 ao 
and w a = 1 ao, respectively, while the system described in 



panel PT2R. which represents the symmetric case w a = Oao, 
is taken as benchmark. The width of the inner well follows 
Eq. (0 with w = 5 ao, shrinking from W lw = 1 ao for 
R = 4a to W iw = 0a for i? = 5a . The width of the 
outer well is assumed constant, W ow = w. The depth of 
the wells is V OW = -15 Hartree and s V™ = -10 Hartree 
for the outer and inner well discontinuity (barrier) rispective. 
Henceforth, the letters A, B, C, D are used to indicate the four 
systems whose potential is depicted in the respective panels of 

Figim 

The entanglement (linear entropy) and its derivatives, 
dL/dR and d 2 L/dR 2 , are plotted in the three panels of 
Fig-El respectively. For R w 4 ao, the many-body wavefunc- 
tion is well-confined inside the inner well. Thus the entangle- 
ment of systems A to D shows there approximately the same 
behaviour, independently of the position of the inner well with 
respect to the symmetry axis of the outer well. The second 
row of panels in Fig Q~2] namely A1,B1,C1 and Dl, shows 
the contour plots of the two-body wavefunction with respect 
to the particles' coordinates x\ and X2 in correspondence of 
the minimum of the linear entropy. With the exception of the 
case D ( Fig.fTZDl ), the contour plots are very similar where 
the wavefunction is significantly different from zero. This, in 
turn, is reflected in a similar position of the minimum of the 
linear entropy. In fact, using the subscripts A, B, C and D to 
indicate the four potentials, we have that the minima of L(R) 
are at R A = 4.52 a , Rb = 4.52 a , Rc = 4.53 a , and 
R D = 4.44 a Q . 

As R increases, the ground-state approaches the migra- 
tion point. Consequently, the particles become progressively 
bounded by very different confinement potentials: symmet- 
ric for A, and increasingly asymmetric from B to D. As the 
ground state wavefunction 'spreads' in the upper well, the 
larger is w a , the more the wavefunction is distorted by a com- 
bination of the pinning effect due to the finite size of the inner 
well and the presence of the left hand side potential barrier of 
the outer well. The third row of panels in Fig Q~2] highlights 
the increasingly marked asymmetry in correspondence of the 
migration points, located at Ro = 4.76 ao, Rc = 4.89 ao, 
R B = 4.92 a and R A = 4.94 a (panels from A2 to D2). 

The most interesting feature is, however, the marked max- 
imum of the entanglement entropy for the three asymmetric 
potentials. In particular, the maximum of Lo{R) is twice 
as large as the maximum value reached in the symmetric 
case. For the symmetric case, the maximum entanglement 
is reached at R = 5 ao. The other maxima are located at 
Rd = 4.76 ao, Rc = 4.90 a and Rb = 4.95 ao, hence 
the maximum of the entropy is reached at a value of R which 
become closer and closer to the migration point as the asym- 
metry increases. For case D, the two values coincide. These 
and the appearance of the entanglement maxima may be un- 
derstood by looking at Fig.[T2"K3 to Fig.[T2"D3. which display 
the contour plot of the wavefunction at the entanglement max- 
ima: (i) the pinning effect localizes most of the wavefunction 
at the location of the inner well, while (ii) the reduced dis- 
tance of outer well walls from the inner well in the asymmet- 
ric potentials further reduces the possibility for the wavefunc- 
tion to significantly expand in the outer well. Due to (i) and 
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FIG. 12: (Color online) Panels B to D: Asymmetric Disappearing inner Well (ADW) potential w a — 1 ao Vs the coordinate x = x\. Panel 
B and C, same as above but for w a = 1.5 ao and w a = 2 ao. Panel A: symmetric case (w a = Oao). Panels Al to Dl: Contour plots of the 
wavefunction against the particles' positions (xi, x-£) at the minimum of L(R). Panels A2 to D2: same as above but in correspondence of the 
'migration' point. Panels A3 to D3: same as above but in correspondence of the maximum of the entropy. 



(ii), as the asymmetry increases, the wavefunction becomes 
more strongly separated into two lobes and more strongly lo- 
calized in two narrow peaks. These two characteristics imply 
that more information could be learned about the second par- 
ticle position once the first is measured and hence means a 
larger entanglement with respect to the symmetric case. In 
particular, for case D, both (i) and (ii) are maximized at the 
migration point, which then corresponds to the maximum of 
the entanglement. 



Using the same argument as in Sec. I Villi it is easy to show 
that the expectation value of the first derivative with respect 
to R of the potential ADW, is discontinuous at R = 5ao- 
It follows that for this value of R, the function L(R) is non- 
differentiable for all four cases, see Fig.[T3b. Finally, we note 
that, for each asymmetric potential, the minimum of the sec- 
ond derivative of the entropy cPL/dR 2 occurs approximately 
in correspondence of the migration point; while the extreme 
point of the ratio (U) / (T) signals the migration point for both 
symmetric and asymmetric geometries of the confinement, see 
Fig.[H 



XI. QPT-LIKE TRANSITION (ASYMMETRIC SYSTEMS) 

Similarly to the symmetric case (see Sec. I VIII) . also for 
asymmetric potentials the minimum of d 2 E/dR 2 and of the 
fidelity F(R, R+SR) is reached at the same value of the driv- 
ing parameter: Rd = 4.72 ao, Rc = 4.86 a , Rb = 4.89 a 
and Ra = 4.91 a , see Fig.fHK and Fig.[T5t. In addition, the 
minimum of dF(Ro ,R)/ dR also occurs at very similar values 
of R, namely R D = 4.72 a , Rc = 4.87 a , R B = 4.90 a 
and Ra = 4.91 a . Notably the minima of d 2 E/dR 2 and 
F(R, R + 5R) become the more pronounced the more asym- 
metric the system is. 

The shoulder structure of dL/dR observed in symmetric 
systems in correspondence to the migration point (and dis- 
cussed in Sec. lVIII ) develops now into a maximum-minimum 
structure (middle panel in Fig. foi l. Here the maximum cor- 
responds to the minimum of d 2 E/dR 2 fairly precisely for 
systems B (maximum of dL/dR at R = 4.90, minimum of 
d 2 E/dR 2 at R = 4.89) and precisely for systems C and D. 
This coincidence of d 2 E/dR 2 and dL/dR extrema has been 
shown to signal QPTs in some systems l5l l43ll . We then ex- 
pect the system at hand to display related characteristics, and 
in particular a rapid evolution between two (almost) orthogo- 
nal states. This is confirmed by Fig. [T5B showing that when 
the driving parameter is swept across the minimum structure 
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FIG. 13: Panel A, B and C: L(R), dL/dR and d 2 L/dR 2 plotted as 
a function of the driving parameter R for the four potentials depicted 
in Fig. [T2] The region of major variations of these quantities are 
detailed in the inset inside the respective panels. The dotted vertical 
lines highlight the migration point for each potential. 
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FIG. 15: Panel A, B and C: d 2 E/dR 2 , F(R ,R) and F(R,R + 
SR) plotted as a function of the driving parameter R for the four 
potentials depicted in Fig. [12] Ro is the value of the parameter at 
which L(R) reaches its minimum for each potential in Fig.fT2] where 
SR — 0.002 ao for all four cases. The inset represent the region 
around the migration point for the panels A and C, while in panel B 
the inset depicts the derivative of F(Ro, R) with respect to R. The 
dotted vertical lines highlight the migration point for each potential. 



tween two almost orthogonal states. 

The minimum of d 2 E/dR 2 and of F(R, R + 8R) becomes 
the more pronounced the more asymmetric the system is, in- 
dicating a more marked QPT-like behavior. We expect then 
a qualitative change of the wavefunction to occur on an even 
smaller parameter range, i.e. when sweeping the driving pa- 
rameter just across the minimum itself. This is indeed the 
case: by looking at the wave-function countour plots for sys- 
tem D, we see that by sweeping the driving parameter in the 
range 4.71 < R < 4.74 the wavefunction evolves from hav- 
ing a single maximum positioned at the inner well, to hav- 
ing additional, marked, lateral maxima. These display the 
increasing importance of Coulomb correlations over confine- 
ment (compare Fig. [T6l panels A and B). 

Interestingly in the asymmetric systems the 'migration' 
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FIG. 16: Panel A, B, C: Asymmetric Disappearing inner Well 
(ADW) case D. Contour plots of the wavefunction against the par- 
ticles' positions (xi, £2) at R — 4.72 ao (minimum of d 2 E/dR 2 
and F(R, R + SR)) R = 4.74 a and R = 4.80 a respectively 



point marks a secondary structure (a shoulder) on the right of 
the minimum of F(R, R + 6R) and d 2 E/dR 2 . In this region 
the wavefunction maximum located above the inner well dis- 
appears, and, due to Coulomb correlations, a valley develops 
between the wavefunction lobes (Fig.[T6C). 

The very narrow parameter region between the minimum 
of F(R,R + 5R) and d 2 E/dR 2 and the 'migration' point 
represents the crossover from confinement to Coulomb corre- 
lations as the leading term in shaping the wave-function char- 
acteristics. In this region the entanglement is extremely sensi- 
tive to variations of the external parameter, going from its al- 
most minimum to its maximum value (seeFig.[T3V This very 
strong sensitivity could be exploited to envisage entanglement 
'switches'. 



XII. CONCLUSIONS 

We studied a set of systems of two electrons confined by 
a potential which evolves from a core-shell to a single well 
potential. For symmetric confining potentials (DIW system) 
as the driving parameter R increases, the inner well width is 
reduced while the width of the outer well increases. The sys- 
tem ground state is initially bounded to the inner well with 
a wavefunction in a quasi-product state (entanglement mini- 
mum and Coulomb energy maximum). In the region around 



the point where the two electrons' ground-state migrates into 
the outer well, the energy and the entanglement entropies are 
continuous and differentiable. However, here various QPT 
markers display a behavior similar to the one of a second- 
order QPT: in particular the second derivative of the ground 
state energy, as well as the fidelity F(R, R + 6R), displays 
a minimum, and the first derivative of the ground state entan- 
glement presents a structure around this minimum, we then 
associate this parameter region to a 'QPT-like' transition. The 
wavefunction indeed undergoes here a drastic change, quickly 
evolving from a quasi-product towards a triplet-like state. We 
showed that a similar QPT-like transition characterizes also 
the evolution from a core-shell to a double-well potential de- 
scribed in Ref . [ij. 

We compared the results for the DIW system with a bench- 
mark system where only the outer well is present, and noticed 
that a very narrow inner well has a strong pinning effect on 
the entanglement (and fidelity) of the system. In particular for 
this symmetric system, the entanglement may be reduced by 
half by the presence of an inner well, even when the ground 
state is already bounded to the outer well. 

We also consider systems (ADW) with the position of the 
core well being asymmetric with respect to the shell. We show 
that, in the 'migration' region, their sensitivity to small varia- 
tion of the driving parameter is even greater than for the sym- 
metric case, with the entanglement displaying a sharp maxi- 
mum at the 'migration' point and a value as high at twice the 
maximum value of a system without the inner well and about 
four times the value of a corresponding DIW system. This 
property might be exploited to create entanglement switches 
by moving the position of the inner well with respect to the 
outer well confining potential. 

The maximum in the entanglement of ADW systems de- 
rives from a combination of the 'pinning effect' by the inner 
core and of the asymmetric confinement (in respect to the pin- 
ning center) provided by the outer shell. This interplay results 
in the high probability of the second particle to be found in 
a more sharply defined small region of space, once the po- 
sition of the first is measured, hence the increase in entan- 
glement. In this region the analysis of QPT markers, and in 
particular the coincidence of the minimum of the energy sec- 
ond order derivative with the maximum of the entanglement 
entropy first order derivative, confirms a QPT-like behavior 
for the system wavefunction. This transition occurs in the re- 
gion where, as the wavefunction expands from the inner to the 
outer well, Coulomb correlations starts to substantially mod- 
ify its shape, and is due to the interplay between confinement 
by the core well, Coulomb interactions, and confinement by 
the outer well. 

We demonstrated that a potential which is characterized by 
a driving parameter R and has a finite discontinuity, even just 
at a single point (x, R), may induce a discontinuity in the 
derivative of the ground-state wavefunction at R for any (or 
almost any) value of x. This in turn induces a non differentia- 
bility at R, with respect to R, of any functional of the ground 
state and in particular - and in contrast with QPT signatures 
- of the same order derivatives of energy and entanglement. 
This is the case for the DIW potential described in this work: 
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here the non-analytic behavior of ground state energy and en- 
tanglement, which is reminiscent of the one encountered in 
QPT, derives instead in a non-trivial way from the finite dis- 
continuity at a single point of the confining potential. To un- 
derline the peculiarity of this connection we also presented 
a counter-example in which similar discontinuities in the po- 
tential do not transfer to other quantities. This is the case of 
the rectangular-like limit of the potential considered in Ref.Ql 
where exact cancellations occur, and hence the ground state 
wavefunction and all related quantities, such as the entangle- 
ment and energy, remain differentiable at any R. 

We presented a detailed analysis of the particle density fi- 
delity, and showed that, for the ground state, this quantity pro- 
vides similar information to the wavefunction fidelity, but may 
be calculated directly from the more accessible particle den- 
sity. In particular we demonstrate that, for ground-states with 
nodeless spatial wavefunctions, the particle density fidelity is 
zero if and only if the wave function fidelity is zero. 

Our results suggest that the entanglement of two particles 
confined in a controllable core-shell structure would display 
QPT-like characteristics as the ground states 'migrates' from 
the inner well to the outer shell, including a well defined min- 
imum of both the wavefunction and particle density fidelities. 
In the corresponding narrow driving parameter range our cal- 
culations demonstrate a high sensitivity of the entanglement 
to variations of the core well size and/or position. 

Finally our results show that the particle density fidelity 
may complement or perhaps even replace the traditional wave- 
function fidelity as a diagnostic tool of the ground state of 
those systems for which the particle density may be experi- 
mentally accessed. 



Appendix A: Effect of the finite discontinuity of the confining 
potential 

1. Core-shell to single well potential (DIW potential) 

We note from Eq. (0 and Eq. (0 that Vdiw(x; R) un- 
dergoes linear transformations (contraction of the inner well, 
dilatation of the upper well) which are continuous every- 
where in R with the exception of the point (x = 0; R = 
w). Treating the potentials Vdiw(x;R) as a family of 
square integrable functions (this condition is verified for 
R < oo), the usual distance can be defined and in partic- 
ular we have that D{Vdiw(x; w — SR) 7 Vdiw(x;w)) = 

\J J \ Vdiw(x, w — 6R) — Vdiw(x, w)\ 2 dx goes to zero 
continuously as SR —> 0. This means that the disappearance 
of the inner well at R = w marks the transition between two 
potentials which differ over a set of measure zero. 

We now consider the expectation value of the derivative of 
our potential, (dVDiw(x;R)/dR), which is directly related 
to the derivative of the ground-state energy, see Eq. i22i . 

Of the two components of Vdiw, see Eq. (0, only the 
derivative of V™ will contribute to a possible discontinuity 
at (x = 0; R = w), as the outer well is simply widening as R 
increases. Its derivative is given by 



dV iw (x;R) _ 
dR ~ 
V \ J w-R\ n ( w~R 



w-R\ r / w-R 
x H — }S \-x + — - — 



(Al) 



We can calculate (dV tw (x; R) / dR) with the help of the prop- 
erty 

n(x; R)6(±x + xq)S(^x + xo)dx = c XQ n(±xo; R) 

-oo 

(A2) 

where c Xo = 1 for xq > 0, c X(l = for xq < 0. By using the 
symmetry n(x; R) = n(—x; R) and considering R < u>, we 
find 



t dV iw {x-R) 
" dR 



V 



n(W iw /2;R). 



R<w 



On the other hand, for R > id, we obtain 

8V iw (xf,R) 



dR 



= 0. 



(A3) 



(A4) 



R>w 



For R — > w 1 , the above equations give 



,dV iw (x;R) 
' dR 



^V^jx-R) 
' dR 



—n(Q-w). 
(A5) 

Here n(0; w) is different from zero, and hence the expecta- 
tion value of the derivative of the potential with respect to R 
is discontinuous at R = w. 

As a consequence, the derivative of the wavefunction with 
respect to R is discontinuous at (x; R = w) for all or al- 
most all x, and this explains the discontinuities observed in 
F(R,R + SR) and F n {R,R + SR) and in the derivatives 
in respect to R of entanglement, energies, F(Rq,R), and 
F n (R 07 R). 



2. Core-shell to double well potential (DW potential) 

Similarly to the DIW potential, the potential in Eq. (0 dis- 
plays a finite discontinuity in respect to R at (x = 0; R = w). 

For R < w, Eq. (0 is similar to the DIW potential except 
that in Eq. (0 the outer well decreases as the inner well does. 
For R> w, Eq. (0 describes two separated wells which move 
further apart and decrease in width as R increases. 

By repeating calculations similar to the ones done for 
Vdiw we obtain 
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{ ) =-—[n(0;w)+n{w;w)],(A6) 

an R ^ w ± 1 

which shows that no discontinuity appears in the limit R — > 
w in the hypothesis that the wavefunction, and hence 
n(x,R) is a continuous function of R: in this case, even if 
dVow /dR does contain delta function-type discontinuities, 
there is no discontinuity in (d Vdw I dR) and hence the deriva- 
tive of the wavefunction is continuous as well as the deriva- 



tives of the other functions discussed in this paper. This pic- 
ture is confirmed by Fig. [10] where calculations performed di- 
rectly with the square-well potential in Eq. © show a steep 
gradient of the entanglement at R = w = 5 eto, but not a 
discontinuity in its derivative. 

We remark that this cancellation of the discontinuity is ac- 
cidental and is due to the fact that the depth of the inner well 
is the same as the height of the barrier between the double- 
well structure. The discontinuity reappears if this symmetry 
is lifted. 
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